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Using neutron powder diffraction and Monte-Carlo simulations we show that a spin-liquid regime 
emerges at all compositions in the diamond-lattice antiferromagnets Co(Ali_a;Coi)204. This 
spin-liquid regime induced by frustration due to the second-neighbour exchange coupling J2, is 
gradually superseded by antiferromagnetic coUinear long-range order (k=0) at low temperatures. 
Upon substitution of Al^^ by Co"^^ in the octahedral B-site the temperature range occupied 
by the spin-liquid regime narrows and Tjv increases. To explain the experimental observations 
we considered magnetic anisotropy D or third-neighbour exchange coupling J3 as degeneracy- 
breaking perturbations. We conclude that Co(Ali_a;Co3;)204 is below the theoretical critical point 
J2/J1 — I/8, and that magnetic anisotropy assists in selecting a coUinear long-range ordered ground 
state, which becomes more stable with increasing x due to a higher efficiency of 0-Co'^"'"-0 as an 
interaction path compared to 0-Al'^"'"-0. 
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I. INTRODUCTION 



Magnetic systems with frustration induced by com- 
peting exchange interactions quite often manifest 
unconventional ground states, the most intriguing of 
which are spin-liquids^. Recently one such exotic state, 
a 'spiral spin-liquid', was uncovered theoretically in a 
classical treatment of diamond-lattice Heisenberg anti- 
ferromagnets (AFM) by Bergman et alA who showed 
that competition between nearest and next-nearest 
neighbor exchange couplings Ji and J2 creates - for 
•h/Ji > 1/8 - a highly degenerate ground state con- 
sisting of a set of coplanar spirals, whose propagation 
vectors form a continuous surface in momentum space. 
The frustration results in a rich phase diagram as a 
function of the ratio J2I J\- The degeneracy of these 
ground states can be lifted by thermal^ or quantum^ 
fluctuations leading to an 'order-by-disorder' phase tran- 
sition from a spiral spin-liquid regime to an ordered state. 

Among the diamond-lattice AFM, compounds with 
the spinel structure recently attracted much attention^^— . 
In particular, Co-Al oxides were considered as promising 
candidates for study of 'order- by-disorder' physics^. In 
these compounds of general stoichiometry AB2O4 the 
tetrahedral A-sites are occupied by high-spin (S=3/2) 
magnetic Co^+ ions which form a diamond lattice 
consisting of two interpenetrating face-centered cubic 
sublattices coupled antiferromagnetically. The octahe- 
dral B-sites can be filled either by nonmagnetic AP+ 
ions and/or by low-spin (S=0) nonmagnetic Co'^"*" ions. 



Existing experiments on Co-Al oxide spinels do not 
provide a clear picture. An early neutron-diffraction 
study^ on C03O4 showed that the magnetic moments of 
the Co^+ ions located at the tetrahedral sites form a sim- 
ple coUinear AFM below the Neel temperature T7v=40 K. 
This picture has been questioned by a recent /j,SR studj*iS 
which found two frequency components near Tat suggest- 
ing incommensurate magnetic order. Experimental ob- 
servations on C0AI2O4 are also contradictory. A powder 
neutron diffraction study of Rotb^i suggested long-range 
AFM order below 4 K, while Krimmel et alw^ detected 
a spin liquid (or glassy-like) ground state. Electron spin 
resonance (ESR), magnetization and specific heat mea- 
surements identified the ground state of C0AI2O4 as spin- 
glass-like with a high frustration parameter ITcwI/Tat 
of 22 [Reflill and 10 [ReflU], respectively. These ex- 
perimental results are consistent with the calculations 
of Bergman et alA which placed C0AI2O4 in the region 
of J2/J1 ~ 1/8 where the spiral surface begins to de- 
velop. To clarify the situation, it is essential to un- 
derstand whether the ground state of C03O4 is a sim- 
ple coUinear antiferromagnet and if C0AI2O4 is a spiral 
spin-liquid fluctuating among degenerate spirals. It is 
also important to understand why substitution within 
the nonmagnetic B-site changes the magnetic properties 
so drastically. Tristan et alr^ attempted to answer these 
questions based on bulk macroscopic measurements of 
Co(Ali_a;Coa;)204 system. They proposed that the spin- 
liquid state is realized in C0AI2O4 (x=0), while with in- 
creasing Co substitution x the second neighbor coupling 
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J2 decreases and collinear AFM long-range order devel- 
ops. 

In this work we address these important issues and refine 
the evolution of the magnetic states in Co(Ali_a; 002^)204 
spinels by means of neutron powder diffraction sup- 
ported by Monte-Carlo simulations^^. Our analysis re- 
veals that for x=0 the system is close to the critical 
point J2/</i=l/8 in the phase diagram while frustration 
is weaker for x>0. We argue that magnetic anisotropy 
assists in the selection of a collinear long-range ordered 
ground state. The difference in Tjv and in the extent 
of the spin-liquid regime apparently originates from the 
decrease of the J2 / Ji ratio with increasing x. 
Polycrystalline Co(Ali_2.Coa;)204 samples with x=0, 
0.35, 0.75 and 1 have been prepared as reported in 
Ref. [13I . Thorough structural characterizations by high- 
resolution x-ray synchrotron diffraction (on the Materials 
Science beamline at Swiss Lig ht Source, A=0.41414 A) 
and neutron diffraction (on the HRPT diffractometer at 
the Swiss Neutron Spallation Source SINQ, A=1.1545 A) 
indicated the absence of inversion in all samples except- 
ing C0AI2O4. For the latter compound the inversion is 
17% and this might infiuence T^. Analyzing the lattice 
constants and x-ray peak profiles we concluded that the 
Co'^+ ions are homogeneously and randomly distributed 
over the B-sites. This follows because (i) the lattice 
constants obey Vegard's law as the function of x and 
(ii) there is no peak asymmetry, which would have been 
present in the high-resolution data if any significant de- 
gree of inhomogeneity was present. 
Medium-resolution neutron powder diffraction patterns 
for four compositions have been collected in the tem- 
perature range 1.5 K - 150 K on the DMC instrument 
at SINQ with a neutron wavelength A=2.4526 A. For 
all compositions, broad bumps resembling a liquid-like 
structure factor and reflecting short-range magnetic or- 
der start to develop below the Curie- Weiss temperature 
^cvf|=110 K— ; they narrow and shift to higher sm9/X 
with cooling (Fig. [T]). Approaching Tjv the diffuse scat- 
tering localizes near the <111> and <200> positions. 
Below Tjv the liquid-like features remain but gradually 
loose spectral weight as magnetic Bragg peaks due to 
long-range order develop. The spectral weight of the dif- 
fuse scattering component continuously increases from 
C03O4 to C0AI2O4 (Fig. [5]) as the frustration parameter 
[TcwI/Tat grows. We can quantify it by the area of a 
Lorentzian fitted to the first diffuse bump. As shown in 
the inset of Fig.[2]this area is largest roughly at T/Tjv ~1 
for all compositions, but for x=0 and x=0.35 diffuse scat- 
tering develops far above and remains significantly below 
this value. For x=l the T/T^r interval revealing diffuse 
scattering is narrower but still much to large to be inter- 
preted as classical critical scattering of a 3D long-range 
ordered magnet. 

Regarding the long-range order, all samples show mag- 
netic Bragg peaks at low temperatures. The ordering 
temperature Tjv and the static ordered magnetic mo- 
ment decrease from C03O4 to C0AI2O4 (see Table |T]and 
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FIG. 1: (Color online) DMC patterns of x=0 and x=l sam- 
ples at selected temperatures. Reasonable statistics on diffuse 
scattering has been obtained by measuring for about 10 h at 
one temperature. Inset: Temperature dependence of intensity 
of the <200> magnetic reflection normalized to the nuclear 
<220> intensity. 



TABLE I: The Neel temperature, Tjv, and the ordered mag- 
netic moment, M, determined at 1.6 K from the DMC pat- 
terns. 



X Tiv [K] M [fiB] 

1.0 29 3.53(3) 

0.75 16.5 2.59(3) 

0.35 9 1.31(4) 

0.0 5 0.25(7) 



inset of Fig.[T]). It should be noted that the x=0 diffrac- 
tion pattern does not correspond to a conventional long- 
range ordered state: diffuse scattering clearly dominates 
and the <200> peak is so broad and weak that the or- 
dered moment cannot be determined with high accuracy. 
Nevertheless, the magnetic Bragg pattern is the same for 
all compositions and is consistent with the collinear two- 
sublattice model proposed by Roth^. 

To model the observed magnetic diffuse scattering 
we used the quasistatic approximationi^, which assumes 
that the observed diffuse scattering can be attributed 
to static correlations. This approach is justified since 
the bandwidth of magnetic excitations is below 6 meV 
whereas all scattering with energies below 14.7 meV is 
summed in our diffraction experiment. We performed 
a Monte-Carlo search for the ground state of a clus- 
ter of 1047 atoms. The energy of the classical Heisen- 
berg Hamiltonian for spins S interacting with first- and 
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FIG. 2: (Color online) Maximal diffuse magnetic scattering 
in the differential T-150 K neutron powder patterns. The T- 
value is equal to 6, 15, 19, 45 K for the compositions x=0, 0.35, 
0.75 and 1, respectively. Inset: Evolution of diffuse scattering 
quantified as an area of Lorentzian fitted to the first bump. 

second- neighbour antiferromagnetic couplings Ji , J2 > 
if = Ji ^ S, • S, + J2 ^» • 

<ij> <<ij>> 

was minimized (< ij > means that the sum runs over 
first- neighbour pairs, << ij » on second- neighbours 
pairs) . If we add either a third- neighbour interaction or 
an anisotropy contribution, the Hamiltonian will be com- 
pleted with an additional term AH, given, respectively, 
as 

AH = J3 • (third-neighbour); (2) 

<<<y>>> 

AH = DYi^f^f (anisotropy). (3) 

i 

where <<< ij »> refers to third- neighbours pairs and 
u is the anisotropy direction, that we took as (111). D 
is the magnitude of the anisotropy term. The moments 
were kept equal and constant in magnitude; their direc- 
tion was changed at random, one at the time, and to ob- 
tain the ground state only energy-decreasing moves were 
accepted in the final stage. The stopping criterion for 
the Monte-Carlo was that the last 1000 accepted config- 
urations had an energy spread less than 10~^ Ji. Runs 
have been repeated for the whole range 0< J2/J1 <1. 
We calculated the static spin-pair correlation functions 



and diffraction patterns for the Monte-Carlo ground state 
(MCGS) and the reference AFM clusters for each J2/ Ji- 
The correlation function is given by 

CF{d)='^-YS,-S,5{\n-r,\-d) (4) 

ij 

where d is the distance between spins at positions 
and rj. We found that for the range 0< <l/8 
the coUinear AFM is the ground state, in agreement 
with Ref. i. Remarkably, the MCGS was always config- 
urationally different from a coUinear AFM state even if 
the energy difference was negligible (order of 1 mK or 
less). Analytical calculations (Fig. [3] inset), similar to 
Ref . 13, support this result. In fact, the energy minimum 
corresponding to the ground state for J2/J1 <l/8 is 
extremely flat around the g=0 point of the first Brillouin 
zone, and therefore, very many states are in the range of 
thermal excitation at any accessible temperature. 
The MCGS starts closely related to the coUinear 
AFM state, but progressively departs from it with 
increasing Jij J\- This is clearly seen in the decay of 
spin correlations, which can be quantified by the ratio 
CFmcgsICFaf- In Fig.[2]we plot \og\CFMCGslCFAF\ 
versus squared distance between spins. Near 
J2I Ji=^/i the decay changes from predominantly Gaus- 
sian e"^/^'^'^/'") (a straight line in this choice of axes) to 
a mixed Gaussian-exponential form e~'^/^ ■ e~'^l'^''y'^l'^^ 
(a curve bending downwards for 200< 0?^ <400 A^). 
This behaviour can be attributed to the incipient spiral 
surface. In the calculated diffraction patterns (Fig. SI 
top) these changes are evidenced as progressive broad- 
ening and distortion of Bragg peaks. Pure Gaussian 
correlation decay yields a Gaussian peak broadening, 
whereas gaussian-exponential decay yields Voigtian 
point spread (that is, with a sizable Lorentzian convo- 
lution component). For J2/J1 >0.5 the correlations are 
completely changed and features of the coUinear AFM 
state fade out. 

We now compare the experimental and calculated 
diffuse scattering. The patterns shown in Fig. [2] were 
obtained as differences between data measured at T= 6, 
15, 19 and 45 K (for x=0, 0.35, 0.75 and 1, respectively) 
and 150 K. In the calculated patterns in order to obtain 
the ordered AFM state at T < T/v we needed to perturb 
the simple Heisenberg Hamiltonian so that the AFM 
minimum would be deeper and the AFM configration 
more stable. For this we considered separately two 
possible degeneracy-breaking perturbations: magnetic 
anisotropy and finite J3 coupling. To retain a simple 
model we described anisotropy as a single-ion parameter 
D. The MC procedure was now changed, allowing for 
system equilibrations at different temperatures. All 
moves that would decrease the energy or that would 
increase it with probability oc e~'^^/'^ were accepted^. 
When the average and Uuctuations of the energy of a 
large number of the last accepted states were sufficiently 
stable the temperature was changed. Fixing Ji=l as 
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FIG. 3: (Color online) Log of ratio of spin-pair correlation 
functions \CFi/CFAFM\ versus for selected J2/J1 ratios. 
Inset: The minimal energy curve for the colhnear antifer- 
romagnet (red) and for the best Monte-Carlo ground state 
(black). 



a convenient energy scale, J2, J3, D or T were varied, 
each in several steps. 
Equally good fits were obtained with J3 and magnetic 
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FIG. 4: (Color online) Top: Calculated powder diffraction 
patterns for several selected J2/J1 ratios at T = K. Bottom: 
Observed (symbols) and calculated (line) powder diffraction 
patterns for x=0 and x=l corresponding to the best model 
parameters. The calculated patterns are corrected for the 
magnetic form factor of Co^+ and polynomial background. 



anisotropy D as perturbations, signifying that the 
available diffraction data are not sufficient to determine 
which of these terms stabilizes the ordered ground 
state. However, we give preference to the magnetic 
anisotropy following the theoretical study of Ref. [2l|. 
The best fits presented in Fig. |4] (bottom) suggest 
the ratio J2/Ji=0.125 for C0AI2O4 and J2/Ji=0.05 
for C03O4. The change of the exchange energy with 
substitution apparently originates from the peculiarities 
of the electronic structure. Band structure analysis^^ 
shows that near the Fermi level in C03O4 there are 
Co'^+ d and oxygen p states, while in C0AI2O4 the Al 
p states are absent from the Fermi level and the weight 
of O p is diminished. This implies that the interaction 
path 0-Co'^+-0 is more effective and the corresponding 
exchange integrals are larger in C03O4. 
Our findings give a natural explanation of the fiSK 
rcsultsi^ on C03O4 as due to spin- liquid physics rather 
than incommensurate magnetic order. Also, based on 
the position of C0AI2O4 in the J2/J1 phase diagram, we 
suggest that even in an ideal sample, with no inversion or 
other perturbation, the ground state would be a collincar 
AFM, though one might need very low temperatures to 
reach it. As the ground state is the colhnear AFM the 
general term 'spin-liquid' and not 'chiral spin-liquid' is 
appropriate for the high-temperature regime above Tjv 
in the title system. 

A single crystal inelastic neutron scattering experiment 
would allow to extract the absolute values of Ji , J2 
and to validate our conclusions. We remark that the 
elaborated approach to fit measured diffuse magnetic 
neutron scattering to Monte- Carlo simulations can be 
easily adapted to other frustrated systems and would be 
useful in justification of an anticipated Hamiltonian. 
In summary we studied the evolution of magnetic 
states in Co(Ali_a;Coa;)204 polycrystalline samples with 
temperature and substitution in the B-site. We observed 
short-range and long-range order for all compositions. 
Employing Monte-Carlo simulations we found that for 
x=0 the system is in the vicinity of the critical point 
J2/ Ji=^/&, where the spiral spin-liquid develops^, 
but stays in the weakly frustrated limit for x>0. We 
also found that replacement in the nonmagnetic B-site 
changes the strength of exchange interactions which 
in turn leads to significant differences in the ordering 
temperatures and in the extent of the spin-liquid regime. 
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